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The GECAM series of satellites utilizes LaBr3(Ce), LaBr3(Ce,Sr), and Nal(TI) crystals as sensitive materials 
for gamma-ray detectors (GRDs). To investigate the non-linearity in the detection of low-energy gamma rays 
and address the errors in the calibration of the E-C relationship, comprehensive tests and comparative studies 
of the three aforementioned crystals were conducted using Compton electrons, radioactive sources, and mono- 
energetic X-rays. The non-linearity test results of the Compton electrons and X-rays demonstrated substantial 
differences, with all three crystals presenting a higher non-linearity for X/y-rays than for Compton electrons. 
Despite the LaBr3(Ce) and LaBr3(Ce,Sr) crystals having higher absolute light yields, they exhibited a noticeable 
non-linear decrease in the light yield, especially at energies below 400 keV. The Nal(T1) crystal demonstrated 
an "excess" light output in the 6-200 keV range, reaching a maximum "excess" of 9.2% at 30 keV in the X-ray 
testing and up to 15.5% at 14 keV during Compton electron testing, indicating a significant advantage in the 
detection of low-energy gamma rays. Furthermore, we explored the underlying causes of the observed non- 
linearity in these crystals. This study not only elucidates the detector responses of GECAM, but also initiates 
a comprehensive investigation of the non-linearity of domestically produced lanthanum bromide and sodium 
iodide crystals. 
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I. INTRODUCTION 


Recent groundbreaking advancements in various branches 
of astrophysics, such as gravitational waves, fast radio bursts, 
and cosmic rays have been observed, enabling a new "multi- 
messenger, multi-wavelength" era in astronomy [1-5]. These 
discoveries emphasize the importance of efficient detection 
methods to further understand high-energy astronomical phe- 
nomena. Transient gamma-ray sources, including gamma-ray 
bursts and magnetar flares, play a vital role in the study of as- 
tronomy [6-11]. 

The Gravitational wave burst high-energy Electromag- 
netic Counterpart All-sky Monitor (GECAM) series, com- 
prising satellites GECAM-A/B, GECAM-C, and GECAM- 
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D, was developed to monitor various high-energy electro- 
magnetic events, such as gamma-ray bursts and magnetar 
flares [12-16]. These satellites employ gamma-ray detec- 
tors (GRDs) that utilize different scintillating crystals such as 
LaBr3(Ce) and LaBr3(Ce,Sr) for GECAM-A/B, and a com- 
bination of LaBr3(Ce), LaBr3(Ce,Sr), and NaI(T]) to validate 
new detector technologies for GECAM-C. The fourth satel- 
lite, GECAM-D, uses Nal(TI) crystals, and is scheduled to be 
launched in early 2024. The main characteristics of the GRDs 
are listed in Table 1. 


GRDs serve as primary detectors in the GECAM pay- 
load, and GECAM-A/B utilizes an innovative solution that 
employs LaBr3 crystals coupled with silicon photomultiplier 
(SiPM) readout technology (Fig. 1) [14, 15, 17, 18]. LaBr3 
crystals are advanced inorganic scintillators known for their 
high light outputs, excellent energy and timing resolutions, 
good energy linearity, and short light-decay times. The SiPM, 
which replaces the conventional photomultiplier tube (PMT), 
offers advantages such as a simple and compact structure, 
ease of miniaturization, and efficient readout capability. 


For GECAM-C (Fig. 2), the GRDs employ both Nal(T1) 
and LaBr3 crystals, which are coupled to the SiPM readout 
arrays [16, 19]. The Nal(Tl) crystal is a high-performance 
traditional inorganic scintillator with excellent luminescence 


Table 1. The main characteristics of the GRD in GECAM series. 


GRD Parameter GECAM-A/B 
Type LaBr3(Ce); LaBr3(Ce+Sr) 
Quantity 25 
Area 45.36 cm? 
Energy range 8-2000 keV 
Energy resolution < 18% @59.5 keV 
Detection efficiency > 50% @8 keV 
Deadtime 4 us 


Dome cabin 


Electrons cabin 


Fig. 1. The detector layout schematic of the GECAM-A/B pay- 
load. On each satellite, the detectors are designed with a modular 
approach, consisting of 25 Gamma-Ray Detector (GRD) modules 
and 8 charged particle detector (CPD) modules. 


properties that provides a good resolution for both X-rays and 
gamma rays. Inorganic scintillators are widely used as the 
preferred choice for high-energy X/y-ray detectors in space 
owing to their versatility in shaping and sizing, stability, re- 
liability, reasonable cost, inclusion of heavy elements, high 
density, and efficient detection capabilities for X/y-rays. 


(b) Bottom dome 


(a) Top dome 


Fig. 2. Detector layout schematic of the GECAM-C payload. 
GECAM-C is composed of two detector domes: top and bottom. 
Each dome is equipped with different types of GRDs, including Nal- 
based gamma-ray detectors (NGRDs) and LaBr3-based gamma-ray 
detectors (LGRDs). 


The crystals used in the GECAM satellite series were ob- 
tained from the Beijing Glass Research Institute. To optimize 
the performance of these detectors, it is critical to understand 
their energy responses. Consequently, we conducted an in- 
depth study involving X-ray, Compton electron, and gamma- 
ray tests on the LaBr3(Ce), LaBr3(Ce,Sr), and Nal(T1) crys- 
tals used in these satellites[20—22]. Our findings indicate that 


GECAM-C GECAM-D 
Nal(TI); LaBr3(Ce); LaBr3(Ce+Sr) Nal(T)D 
12 5 

45.36 cm? 103.87 cm? 

15-4000 keV 20-1000 keV 
< 18% @59.5 keV < 25% @59.5 keV 
> 75% @15 keV > 60% @20 keV 
4 us 4 us 


the non-linearity of the three crystals varied when exposed to 
distinct excitation sources. The LaBr3(Ce,Sr) crystal exhib- 
ited the strongest linear response to Compton electrons in the 
low-energy range, whereas the Nal(T1) crystal demonstrated 
the best linear response to X-rays. 

Consistent with previous studies regarding the non- 
linearity of iodide crystals [23, 24], domestically produced 
Nal(TI) crystals exhibited a light yield "excess" phenomenon, 
indicating unexpected advantages in the detection of low- 
energy gamma rays. These insights not only contribute to a 
better understanding of the detector response of the GECAM 
series, but also provide invaluable information for evaluating 
the performance of these domestically produced scintillating 
crystals in the low-energy range of 3-400 ke V[25]. Manufac- 
turers can refer to this study to enhance their understanding 
of crystal non-linearity, potentially facilitating the optimiza- 
tion and improvement of the crystal growth processes and 
doping ratios. Furthermore, this study presents the issue of 
non-linearity in crystals for low-energy gamma-ray detection, 
which is highly significant for addressing errors in detector 
calibration related to the energy-channel (E-C) relationship. 
Based on the results of this study, we will further investigate 
the intrinsic resolution of crystals in future studies. 


II. EXPERIMENTAL SETUP AND TEST PROCEDURE 
A. The Wide-Angle Compton Coincidence Technology 


Figure 3 presents the wide-angle Compton coincidence 
(WACC) experimental setup, which primarily comprises a 
radioactive source, a high-purity Germanium (HPGe) de- 
tector, the scintillation detector under examination, and a 
subsequent data acquisition system[24, 26-30]. LaBr3(Ce), 
LaBr3(Ce,Sr), and Nal(T1) cylindrical samples with diame- 
ters of 25.4 mm were selected for this study. Silicone oil was 
used to couple the encapsulated crystals to the PMTs (R6233- 
100, Hamamatsu Photonics, Japan)[31, 32]. According to the 
user manual of the BE2020 planar germanium spectrometer 
manufactured by Canberra, the HPGe crystal had a thickness 
of 20 mm and a volume of 40,000 mm?, allowing for an en- 
ergy detection range of 3 keV—3 MeV[33, 34]. Based on the 
experimental data, the energy resolutions, indicated by the 
full width at half maximum (FWHM), of the HPGe detector 
used in this study were determined to be 1.58 keV (©°Co, 1.33 
MeV) and 1.15 keV (13’Cs, 662 keV). 

The experiment involved placing a 1°’Cs_ radioactive 


source 


i| scintillator | 


Fig. 3. Schematic diagram of experimental setup for obtaining 
Compton electrons with a wide energy range. 


source at a quarter-circle position around the center of the 
crystal, with a distance of 13 cm between them, Gamma pho- 
tons were emitted by the radioactive source via radioactive 
decay and underwent Compton scattering when they struck 
the crystal. Consequently, Compton electrons were gener- 
ated and absorbed by the crystal, whereas certain scattered 
photons escaped from the crystal and were absorbed by the 
nearby HPGe detector. The distance between the tested crys- 
tal and HPGe detector was maintained at approximately 15 
cm[26]. Lead blocks were positioned between the 137Cs 
source and HPGe detector to provide shielding and minimize 
the incidence of primary gamma photons directly irradiating 
the HPGe detector. Coincidence events across a broad energy 
range were obtained by adjusting the position of the radioac- 
tive source and varying the angle between the source, crystal, 
and HPGe detector. 

A desktop waveform acquisition device with a 10 bit 
@ 2 GS/s (interleaved) or 1 GS/s, which was a DT5751 
digitizer[35], was utilized in this experiment to collect sig- 
nals from the crystal and HPGe detectors (Fig. 4). The HPGe 
detector signal operating at +3500 V underwent shaping and 
filtering using an ORTEC 572A amplifier before being sent to 
channel-0 of DT5751. The output signal from the PMT an- 
ode operating at +1300 V was routed to channel-1 of DT5751 
after the photoelectrons were multiplied by the dynodes[36]. 
The signals from the crystal and HPGe detectors successively 
underwent low-threshold discrimination, delayed stretching, 
and a logical coincidence. The resulting coincidence output 
signal was utilized as an external trigger for DT5751, which 
recorded the corresponding coincidence events and generated 
two data files when triggered externally. The secondary par- 
ticles produced by Compton scattering were absorbed by the 
two detectors in a specific order in the time sequence. For a 
"true coincidence event,” the waveform signal corresponding 
to the crystal appeared before the HPGe detector. 

Figure 5 presents the coincidence matrix that represents all 
the collected data of the events. In Fig. 5 (a), the horizontal 
axis represents the energy deposited on the HPGe detector, 
whereas the vertical axis corresponds to the energy deposited 


on the LaBr3(Ce,Sr) crystal. As indicated in the Compton 
scattering formula (Equation 1)[24], as the incident angle (6) 
of the gamma photon increased, the energy of the Compton 
electrons in the LaBr3(Ce,Sr) crystal also increased. 


a (1) 


Mec? 
1+ E,(1—cos@) 


B, = E,-E,= 


where E, is the energy of the gamma ray radiating from the 


source, Ee is the energy of the Compton electron, E, is the 
energy of the scattered photon, 0 is the Compton scattering 
angle, and mec? is the remaining mass energy of the electron. 
Five scattering angles of 0 were chosen during the experiment 
to obtain data over a broad energy range. 

The diagonal points presented in Fig. 5 demonstrate the 
"true coincidence events" of interest. Each point corresponds 
to a specific scattering angle, and the combined deposited en- 
ergies in both the crystal and HPGe detector are constant at 
661.6 keV. The uneven "spread" among the diagonal at dif- 
ferent energy levels is owing to the diverse energy resolutions 
of the crystal for Compton electrons. Furthermore, the non- 
linear response of the crystal determines the “linearity” of the 
diagonal. An analysis of the coincidence matrix enabled the 
extraction of the energy resolution and non-linear response of 
the crystal for Compton electrons. 

Figure 5 presents the horizontal and vertical lines indicat- 
ing the accidental coincidence events simultaneously detected 
by both detectors. The horizontal line represents the finite 
resolution of the crystal and the vertical line represents the 
excellent resolution of the HPGe detector. The other points 
on the graph denote events in which partial energy was de- 
posited in the detector or where detection occurred after scat- 
tering through the surrounding materials. 

The WACC method accurately measures the energy re- 
sponse and resolution of the crystal detector to the Compton 
electrons. Before the Compton experiments, it was neces- 
sary to calibrate the E-C relationship of the HPGe detector, 
which can be obtained from the energy spectra of multiple 
radioactive sources or directly using the vertical lines in the 
coincidence matrix. 

The HPGe detector offers an outstanding energy resolu- 
tion, making it an excellent standard detector. The energy de- 
posited in the crystal was calculated by subtracting the scat- 
tered photon energy in the HPGe detector from the known 
gamma-ray source energy. In actual data processing, the cut 
width of the HPGe energy axis must be determined based on 
the statistics of the Compton scattering event. Within this 
range, the central value is considered to be the energy de- 
posited in the HPGe detector, and Equation 2 is used to cal- 
culate the energy deposited in the crystal. In this study, cut- 
ting was performed along the HPGe energy axis and projected 
onto the crystal axis. 


< Escin >= E— < Eupce >. (2) 


where E, is the known gamma-ray source energy, 
<E y PGe> is the deposited energy in the HPGe detector, and 
<Escin> is the deposited energy in the crystal. 
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Fig. 4. Diagram of the data acquisition system for coincidence events. 


To understand the effect of the cut width or energy 
window width, we measured the energy resolution of the 
LaBr3(Ce,Sr) crystal to 46.6 keV Compton electrons for dif- 
ferent cut widths. The results indicated that the energy resolu- 
tion remained reasonably stable until a cut width of 4 keV was 
reached (Fig. 6). Wider cut widths led to a broadened scatter- 
ing angle range of relevant valid events and an increase in the 
FWHM of the Compton electron spectrum. Subsequently, the 
resolution deteriorated. The energy resolution of the HPGe 
detector was within 1—2 keV, which must be considered when 
determining a reasonable cut width. It is also essential to en- 
sure a sufficient number of events. Therefore, a cutoff width 
of 4 keV was used when the energy deposited in the HPGe 
detector was less than 615 keV. When the deposited energy 
deereasedto, 615—-661.6 keV, a cutoff width of 2 keV was se- 
lected. 


Multiple truncations of the HPGe energy axis were per- 
formed to obtain the crystal spectra for various Compton elec- 
tron energies. Figure 7 illustrates an example of this ap- 
proach, in which a data range of events from 614 to 616 keV 
was considered at an HPGe energy of 615 keV with a cut 
width of 2 keV to produce the Compton electron spectrum 
(Fig. 7 (b)). A Gaussian-shaped, single-energy electron peak 
was visible and fitted with a Gaussian function, returning an 
energy resolution of 15.81 + 0.25% for the 46.6 keV Comp- 
ton electrons in the LaBr3(Ce,Sr) crystal. 


Atoms or molecules are excited when incident particles de- 
posit energy in a crystal, leading to the emission of scintil- 
lation photons with wavelengths similar to those of visible 
light[37]. The light yield, defined as the number of scintil- 
lation photons per unit of energy deposited in the crystal, is 
described by Equation 3. 


ADC 


"= apl E 


(3) 


where S is the light yield of the crystal, ADC repre- 
sents the peak position of the spectrum after deducting the 
baseline, E is the deposited energy in the crystal, and the 
ADCspe = 8.0321 channels denote the single-photoelectron 
response of the Hamamatsu R6233-100 PMT at a high volt- 
age of +1300 V. The response was calibrated using the LED- 
triggered charge method[38, 39]. 


B. Measured by Radioactive Sources 


Radioactive sources of !8°Ba, 137Cs, 241Am, !52Eu, and 
207Bi were employed across a range of y-ray energies from 
30.85 keV to 1063.7 keV to investigate the gamma-ray re- 
sponses. The tested crystal was coupled to a Hamamatsu 
R6233-100 PMT via silicon oil, and the digitizer DT5751 
acquired the signal waveforms in the self-triggering mode. 
ROOT, a data analysis framework conceived by the European 
Organization for Nuclear Research (CERN)[40], was used to 
analyze the experimental data, including baseline subtraction, 
fitting of the full-energy peak, and analysis of the peak posi- 
tion and FWHM. 


C. Single Energy X-ray Measurements Using the Hard X-ray 
Calibration Facility 


We employed two sets of hard X-ray calibration facili- 
ties (HXCF, Fig. 8) established by the National Institute of 
Metrology (NIM) in Beijing Changping, China[41—43] to in- 
vestigate the energy responses of the three aforementioned 
crystals to X-rays in the 8-120 keV range. The HXCF, which 
plays a substantial role in the calibration of gamma-ray de- 
tectors on GECAM, CubeSats, and SVOM satellites[44—46], 
was first built for the high-energy telescope of the HXMT as 
a calibration facility[47] and comprises four primary compo- 
nents: X-ray generator, monochromator, collimator, and stan- 
dard detector. To shield stray light from the X-ray generator, 
the collimator features apertures of various sizes at the en- 
trance and exit. A low-energy HPGe detector (Canberra In- 
dustries) was used as a standard. Before testing, we calibrated 
the HPGe detector for energy linearity, resolution, and detec- 
tion efficiency using various standard radioactive sources[48]. 

The entire set of testing equipment, including the data- 
acquisition system, was placed inside an X-ray testing cham- 
ber (Fig. 9), and remotely controlled for data retrieval in a 
control room. The energy and flux of the X-rays were deter- 
mined using an HPGe detector. The testing procedures are 
shown in Fig. 10. We used GENIE 2000, which is a spec- 
troscopic data acquisition and analysis software, to record the 
spectral data from the HPGe detector. The crystal detector 
was coupled to a PMT (Hamamatsu Model CR160) using sil- 
icone oil. The signals from the crystal detector were collected 
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Fig. 5. Two-dimensional spectrum of coincidence events in both 
HPGe and (a) LaBr3(Ce,Sr), (b) LaBr3(Ce), (c) Nal(Tl) with the 
137s source, that is, the coincidence matrix. 


using a digitizer (DT5751) and analyzed using computer soft- 
ware for the corresponding spectra. 


The range of X-ray testing was between 8-120 keV in this 
study, with fine measurements of the absorption edge of the 
crystal at a step size of 0.1 keV. The performance of the crys- 
tal detector gradually changed as the X-ray energy increased, 
allowing for a reduced number of test energy points. Owing 
to testing at room temperature (22—23 °C), the detector noise 
was slightly higher, limiting the starting test energy points 
to the range of 8-10 keV. For the two LaBrs3 crystals, the 
PMTs coupled to them operated at —800 V, whereas that for 
the Nal(TI) crystal operated at — 1000 V. 
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Fig. 6. Resolution of Compton electrons in LaBr3(Ce,Sr) versus the 
cut width in HPGe for a 1?’Cs source. 


HI. RESULTS AND DISCUSSION 
A. Light Yield Non-linearity to Compton Electrons 


The light yields of the LaBr3(Ce), LaBr3(Ce,Sr), and 
Nal(T1) crystals were normalized to "1" at an energy of 662 
keV. Figure 11 demonstrates the non-linearity of the light out- 
put of the LaBr3(Ce), LaBr3(Ce,Sr), and NalI(TI) crystals to 
the Compton electrons within the energy range 3-400 keV. To 
better quantify the non-linearity of these crystals, we intro- 
duced a metric known as the non-linearity standard deviation 
(NLSD), denoted by Equation 4, where x; represents the rel- 
ative light yield at each energy point. The NLSD values for 
the LaBr3(Ce), LaBr3(Ce,Sr), and Nal(TI) crystals were 0.11, 
0.03, and 0.06, respectively. A larger NLSD value indicates a 
more significant non-linearity of the crystal. 


For both types of LaBr3 crystals, as the energy of the 
Compton electrons decreased, the non-linearity of the light 
yield gradually increased. Within the measured electron en- 
ergy range, the LaBr3(Ce,Sr) crystal exhibited a better lin- 
earity than that of the LaBr3(Ce) crystal, particularly at ener- 
gies below 20 keV. We hypothesized that the doping of Sr?* 
may have improved the internal energy transfer mechanism 
within the LaBr3(Ce,Sr) crystal, enhancing the energy trans- 
fer efficiency in the low-energy region, thereby ameliorating 
non-linearity. Both crystals exhibited a 10% "defect" light 
output at approximately 5 and 20 keV, respectively. The min- 
imum measurable energy point using WACC was 3.1 keV, at 
which the LaBrs(Ce,Sr) crystal exhibited a “defect” of ap- 
proximately 24%, whereas that of the LaBr3(Ce) crystal ex- 
hibited 35%. This experiment validated the simulation results 
presented by Chae, et al.[17], which demonstrated the non- 
linearity of the electrons, while also affirming the accuracy 
and rationality of both the model and experimental work con- 
ducted by the GECAM research team. 


Projected spectrum of LaBr3-Ce,Sr with HPGe gated at 605 keV by irradiation of 137Cs gamma source 
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Fig. 7. Projected spectra of the LaBr3(Ce,Sr) crystal with the HPGe 
energy gated at (a) 605 keV, (b) 615 keV, (c) 625 keV, and (d) 635 
keV, respectively; that is, the (a) 56.6 keV, (b) 46.6 keV, (c) 36.6 keV, 
and (d) 26.6 keV Compton electron energy spectrum. 


In comparison to the two LaBrs crystals, the NaI(T1) crys- 
tal did not exhibit the monotonic "defect" of the lumines- 
cence phenomenon as the energy of the Compton electrons 
decreased. At approximately 14 keV, the Nal(T1) crystal 
reached its maximum light yield, exhibiting an "excess" light 
output of approximately 15.5%. Beyond 14 keV, the light 
yield gradually decreased as the energy increased. Con- 
versely, as the energy decreased to below 14 keV, the light 
yield decreased. The lowest test energy point was 4.1 keV, 
at which the Nal(TI1) crystal demonstrated a luminosity non- 
linearity “defect” of approximately 14%. 


B. The Absolute Light Yield of Crystals 


The three crystals were irradiated using multiple radioac- 
tive sources to obtain the energy spectra of each crystal for 
different sources. The single-photoelectron responses of the 
Hamamatsu R6233-100 PMT used in the measurements were 
calibrated using the LED-triggered charge method at various 
voltages, enabling the calculation of the absolute light yields 
of the three crystals. The absolute light yields and energy 
resolutions of the samples tested at 661.6 keV are listed in 
Table 2. 


Table 2. The absolute light yield and energy resolution of crystals 
for the 661.6 keV gamma rays. 


Crystal Type Size Light Yield(ph/MeV) Energy Resolution(%) 
LaBr3(Ce) 1"x1" 74196 3.00+0.02 
LaBr3(Ce,Sr) 1"x 1" 63939 3.03+0.01 
Nal (T1) 1"x1" 45445 7.18+0.07 
C. Energy Resolution 
Figure 12 illustrates the energy resolution of the 


LaBr3(Ce), LaBr3(Ce,Sr), and NaI (TI) crystals for the Comp- 
ton electrons in the 3—400 keV range. The energy resolution 
of the NaI(TI) crystal was comparable to that of the LaBr3 
crystals at 16-30 keV (Fig. 12). 


The energy resolution of the crystals was expressed us- 
ing the FWHM of the full-energy peak of the X-ray. Fig- 
ure 13 presents the energy resolution of the LaBr3(Ce,Sr), 
LaBr3(Ce), and Nal(TI) crystals for X-rays in the 8-100 keV 
range, as measured by HXCE. The LaBrsz(Ce,Sr) crystal ex- 
hibited the best energy resolution within this energy range. 
At 100 keV, the resolution of the LaBr3(Ce,Sr) crystal was 
8.74 + 0.0681%, whereas those of the LaBr3(Ce) and Nal(T1) 
crystals were 9.41 + 0.0976% and 10.39 + 0.1168%, respec- 
tively. Furthermore, a slight degradation in the energy reso- 
lution, less than 1%, was observed near the binding energy of 
the K-shell electrons. 
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Fig. 8. Hard X-ray Calibration Facility. Both the HPGe and crystal detectors were placed on a displacement platform and maintained on the 
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Fig. 9. Data acquisition system diagram for single-energy X-ray detection. 
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Fig. 10. Test procedure for the non-linear response of the crystals to 
X-rays. 


D. Comparison of the X/y-Ray and Compton Electron 
Responses 


All the data in this study were standardized by setting the 
full-energy peak response of the 662 keV gamma rays from a 
137Cs source as the normalization factor. The non-linearity of 
the LaBr3(Ce,Sr) crystal light yield for the Compton electrons 
and gamma rays in the 3—1000 keV range is shown in Fig. 14. 
Note, the response of the Compton electrons exhibited an ex- 
cellent linearity at approximately 70 keV, with a non-linearity 
of less than 2%. However, a "deficiency" in the light output 
occurred when the energy of the Compton electrons was be- 
low 70 keV, whereas substantial non-linearity was observed 
in the response to gamma rays below approximately 200 keV. 


A more detailed test was conducted on the photon response 
below 120 keV using HXCF. Figure 14 presents the non- 
linear light yield response curve of the LaBrs(Ce,Sr) crys- 
tal to X-rays in the energy range of 8-120 keV. Because 
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Fig. 11. Comparison of the light yield non-linearity of LaBr3(Ce), 
LaBr3(Ce,Sr), and NaI(TI) crystals for the 3-400 keV Compton 
electrons. 
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Fig. 12. Comparison of the energy resolution of the LaBr3(Ce), 
LaBrs(Ce,Sr), and Nal(T1) crystals for the Compton electrons in the 
3—400 keV range. 


the error bars are similar in size to the data point symbols, 
they are not visible in the figure. Ideally, the relative light 
yield should be "1" at all energy points. However, this was 
not the case and varying degrees of light-yield deficiencies 
were observed within the tested energy range. Below 40 keV, 
the LaBr3(Ce,Sr) crystal exhibited substantial non-linearity 
in the relative light yield response to X-rays, with the non- 
linearity exceeding 10%. As the energy decreased, the slope 
of the curve increased, reaching a non-linearity of 36% at 
8 keV. When the X-ray energy exceeded 40 keV, the non- 
linear curve approached the ideal state and the slope became 
milder, indicating that the fluctuation in the number of pho- 
tons generated per unit energy absorbed by the LaBr3(Ce,Sr) 
crystal was small within the energy range of 40-120 keV. 
The LaBr3(Ce,Sr) crystal exhibited absorption edges at 13-15 
keV and 38—40 keV, and a slight reduction in the relative light 
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Fig. 13. Energy resolution of the LaBr3(Ce,Sr), LaBr3(Ce), and 
Nal(TI) crystals varies with the X-ray energy. 
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Fig. 14. Light yield non-linearity of the LaBr3(Ce,Sr) crystal for 
Compton electrons, gamma rays, and X-rays. 


yield was observed within these two energy intervals. 

The NLSD values for testing the LaBr3(Ce,Sr) crystal with 
X-rays and Compton electrons were 0.17 and 0.03, respec- 
tively. The light output of the LaBr3(Ce,Sr) crystal exhibited 
a greater non-linearity in response to the X-rays than to the 
Compton electrons. This can be attributed to the different 
mechanisms by which these particles interact with atoms in 
matter. For X/y-rays ranging from a few keV to several hun- 
dred keV, there are two possible interaction processes with 
the crystal: (1) a direct photoelectric cascade sequence or 
(2) Compton scattering followed by a photoelectric cascade 
sequence. These processes generate several primary elec- 
trons (e.g., Compton electrons and primary photoelectrons) 
and multiple secondary electrons (e.g., Auger electrons and 
secondary photoelectrons), with the final light emission being 
the sum of the contributions from secondary electrons with 
different energies. Notably, these electrons are the products 
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Fig. 15. Light yield non-linearity of the LaBr3(Ce) crystal for Comp- 
ton electrons, gamma rays, and X-rays. 


of the interaction between the incident photons and matter, 
and their energies cannot exceed those of the incident parti- 
cles. Therefore, the light output induced by the photons in 
the LaBr3(Ce,Sr) crystal was always lower than that caused 
by the Compton electrons with equivalent energies. 

We conducted detailed testing of the LaBr3(Ce) crystal us- 
ing the same experimental procedures and data processing 
methods. Figure 15 illustrates the non-linear light-yield re- 
sponse of the LaBr3(Ce) crystal to Compton electrons and 
gamma rays in the 3-1000 keV energy range. The non- 
linear curves tended to be flat, and the results were similar 
for the Compton electrons and gamma rays when the en- 
ergy was above 200 keV; however, significant differences 
were observed below 200 keV. As the energy decreased, the 
LaBr3(Ce) crystal exhibited a lower response to the full- 
energy peak of the gamma rays than to the Compton electrons 
of the same energy. This finding is consistent with the test re- 
sults for the LaBr3(Ce,Sr) crystal, indicating that the manner 
in which the particles interact with the matter directly affects 
the light output of the crystal. For gamma rays in the energy 
range of several hundred kiloelectronvolts, Compton scatter- 
ing is most likely the initial interaction, and most gamma rays 
require multiple interactions for full absorption. The high- 
energy primary and secondary electrons resulting from these 
interactions exhibited good linear responses, reflecting the ex- 
cellent linearity of the response to high-energy gamma rays. 

Figure 15 demonstrates the non-linearity curve of the 
LaBr3(Ce) crystal to X-rays in the 8-120 keV energy range. 
Compared to the LaBrs(Ce,Sr) crystal, this response curve 
deviated from the ideal state more significantly, and nearly all 
the measured energy points exhibited scintillation responses 
below 90%. The light output sharply decreased near the K- 
shell binding energies (13-15 keV and 38-40 keV) of Br and 
La, leading to a greater non-linearity of the LaBr3(Ce) crys- 
tal response curve to the X-rays. The data points below 28 
keV exhibited a non-linearity greater than 20%, and the light 
output at 8 keV was only 58% of the ideal state. 
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Fig. 16. Light yield non-linearity of the Nal(T1) crystal for Compton 
electrons, gamma rays, and X-rays. 


As the X-ray energy decreased, the induced secondary 
electron energies in the crystal decreased, thereby resulting 
in more significant light "defects". The NLSD values ob- 
tained by testing the LaBrs(Ce) crystal with the X-rays and 
Compton electrons were 0.22 and 0.11, respectively. The 
LaBr3(Ce) crystal exhibited a greater non-linearity toward the 
X-rays and Compton electrons compared to the LaBr3(Ce,Sr) 
crystal, particularly at energies below 100 keV. This may be 
attributed to the doping process, indicating that doping with 
Sr?* ions can improve the non-linearity of LaBr3 crystals. 

To better understand the differences in non-linearity among 
the different types of crystals, the Nal(T1) crystal was chosen 
as the third test subject (Fig. 16). Unlike the two LaBr3 
crystals, the NaI(T1) crystal exhibited a pronounced "excess" 
response to the Compton electrons in the 8—80 keV energy 
range, with a non-linearity exceeding 4%. At electron ener- 
gies lower than 6 keV, the crystal displayed slight "defects" 
in the light output, whereas above 80 keV, the curve tended 
to flatten, indicating a good linear response of this sodium 
iodide compound to high-energy electrons. 

Figure 16 also demonstrates the non-linearity of the light 
yield of the Nal (TI) crystal to X-rays in the 8—120 keV en- 
ergy range. Compared with the response to the Compton elec- 
trons, the X-ray test results exhibited a similar trend, with 
NLSD values of 0.06 for both. However, there were differ- 
ences in the slopes of the curves. Direct photoelectric inter- 
actions with matter are most likely to occur for photons in 
the tens-of-keV range. Assuming this photoelectric absorp- 
tion occurs with iodine K-shell electrons (probability of 83% 
when the photon energy is greater than 33.17 keV), the re- 
sulting photoelectrons have energies within the range, with a 
substantial "excess" light output. The total light emission in- 
duced by all secondary electrons generated from the photons 
exceeded that caused by Compton electrons with equivalent 
energies. Therefore, when the energy was within the 40-70 
keV range, the Nal(T1) crystal exhibited a higher relative light 
output to the X-rays, producing a greater number of photons 


per unit of X-ray-deposited energy compared to the case of 
the Compton electron incidence. 

The response of the Nal (TI) crystal to the X-rays was sim- 
ilar to that of the Compton electrons at approximately 33 keV, 
is related to the binding energy (33.17 keV) of the iodine K- 
shell electrons, as photons with energies lower than this can- 
not excite K-shell electrons from the iodine atoms. Nearly all 
of the photon energy was transferred to electrons, and only 
a small fraction of low-energy photons interacted with an io- 
dine L-shell electron (with a binding energy of 5.19 keV) to 
produce lower-energy X-rays through the photoelectric effect. 

Within the measured X-ray energy range, the Nal(TI) crys- 
tal exhibited varying degrees of "excess" light output, which 
can also be explained by the sequence of the photoelectric ef- 
fect cascade. The low-energy electron response demonstrated 
an "excess" and reached its maximum value at 14 keV, as 
shown in Fig. 16. Therefore, when photons undergo a se- 
ries of interactions to produce multiple low-energy secondary 
electrons, a "burst" phenomenon occurs in the light output, 
which also explains why the photon response reached a max- 
imum value at approximately 30 keV instead of 14 keV. As 
the incident photon energy increased, the light output grad- 
ually decreased but remained above 100%, which was ow- 
ing to the more complex distribution of secondary electron 
energies, resulting in numerous secondary electrons with en- 
ergies lower than 6 keV. The electron response below 6 keV 
exhibited a "deficient" luminous response, which formed a so- 
called "compensation" effect with the "excess" phenomenon 
observed in the tens of keV range. 


IV. CONCLUSION 


We employed the WACC technique and HXCF/radioac- 
tive sources to compare the energy responses of domestically 
produced LaBr3(Ce), LaBr3(Ce, Sr), and Nal(T1) crystals 
to Compton electrons and X/y-rays. The NLSD values ob- 
tained through the X-ray testing of the LaBr3(Ce), LaBr3(Ce, 
Sr), and Nal(TI1) crystals were 0.22, 0.17, and 0.06, respec- 
tively. In comparison, the Compton electron testing resulted 
in NLSD values of 0.11, 0.03, and 0.06 for the same crys- 
tals. The non-linear curves of these domestic crystals exhib- 


[1] B. Abbott, R. Abbott, T. Abbott, et al., Observation of gravita- 
tional waves from a binary black hole merger. Phys. Rev. Lett. 
116, 061102 (2016). Doi: 10.1103/PhysRevLett.116.061102 
B. Abbott, R. Abbott, T. Abbott, et al., Gravitational waves and 
gamma-rays from a binary neutron star merger: Gw170817 
and grb 170817a. Astrophys. J. Lett. 848, L13 (2017). Doi: 
10.3847/2041-8213/aa920c 

B. Abbott, R. Abbott, T. Abbott, et al., Gw170817: observa- 
tion of gravitational waves from a binary neutron star inspi- 
ral. Phys. Rev. Lett. 119, 161101 (2017). Doi: 10.1103/Phys- 
RevLett.119.161101 

[4] B. Abbott, R. Abbott, T. Abbott, et al., Multi-messenger obser- 


[2 


fae 


[3 


= 


10 


ited different slopes (Fig. 14, Fig. 15 and Fig. 16), indicating 
varying degrees of non-linearity at low energies. Based on the 
experimental results, the non-linearity of the three crystals to 
the X/y-rays exceeded that of the Compton electrons, which 
can be attributed to the distinct interaction mechanisms be- 
tween the incident particles and material. 

The NLSD values for LaBr3(Ce) were 1.29 times higher 
for X-rays and 3.67 times higher for Compton electrons com- 
pared to LaBrs(Ce, Sr), indicating that the LaBr3(Ce, Sr) 
crystal exhibited a better linearity and that doping with Sr?* 
ions can improve non-linearity. However, the absolute light 
yield of the LaBr3(Ce, Sr) crystal was slightly lower than 
that of LaBr3(Ce) (Table 2), potentially owing to the need 
for further optimization of the growth process and doping ra- 
tio by domestic manufacturers. The energy resolution of our 
LaBr3(Ce, Sr) crystal was inferior to that reported by its for- 
eign counterparts[49]. This discrepancy may have occurred 
owing to the inherent performance variations among the dif- 
ferent crystals, differences in measurement methods when 
coupled with PMT, or distinctions in the growth processes and 
raw materials between the Chinese and Saint-Gobain crystals. 

The Nal(T1) crystal exhibited an "excess" light output of up 
to 9.2% when tested with X-rays and 15.5% when tested with 
Compton electrons. This "excess" light output distinguished 
the Nal(T]) crystal for having the advantage of detecting low- 
energy X/y-rays. The calibration and in-orbit performance of 
GECAM-C validated that the NalI(Tl) crystals exceeded ex- 
pectations [16, 19]. Energy resolution is not the primary con- 
cern in gamma-ray burst detection. Although Nal(TI) crystals 
may not match the energy resolution and absolute light yield 
of LaBr crystals, the test results demonstrated their satisfac- 
tory performance in the 10-1000 keV energy range. Further- 
more, Nal(TI) crystals can be manufactured in large sizes and 
are cost-effective. Consequently, GECAM-D utilizes Nal(Tl) 
crystals as sensitive detection materials. 

We conducted a study regarding the light yield and non- 
linearity of three crystals produced by the Beijing Glass Re- 
search Institute, and found that different calibration standards 
are required for the detection of gamma-rays and electrons. 
While the current GECAM satellite GRDs lack electron- 
gamma discrimination capabilities, the non-linearity results 
of the Compton electrons may be applied to future correc- 
tions for electron detection. 


vations of a binary neutron star merger. Astrophys. J. Lett 848, 
L12 (2017). Doi: 10.3847/2041-8213/aa91c9 

[5] T. Li, S. Xiong, S. Zhang, et al., Insight-hxmt observations 
of the first binary neutron star merger gw170817. Sci. China.- 
Phys. Mech. Astr. 61, 1-8 (2018). Doi: 10.1007/s11433-017- 
9107-5 

[6] V. Connaughton, M. Briggs, A. Goldstein, et al., Localiza- 
tion of gamma-ray bursts using the fermi gamma-ray burst 
monitor. Astrophys. J. Suppl. S. 216, 791-804 (2015). Doi: 
10.1088/0067-0049/2 16/2/32 

[7] M. Stanbro, M. Briggs, O. Roberts, et al., A fermi gamma- 
ray burst monitor event observed as a terrestrial gamma-ray 


[8] 


[9] 


[10] 


[11] 


[12] 


[13] 


[14] 


[15] 


[16] 


[17] 


[18] 


[19] 


[20] 


[21] 


[22] 


[23] 


flash and terrestrial electron beam. J. Geophys. Res.-Space. 
124, 10580-10591 (2019). Doi: 10.1029/2019JA026749 

X. S., X. S.L., C. C., et al., Energetic transients joint analy- 
sis system for multi-instrument (etjasmin) for gecam — i. posi- 
tional, temporal, and spectral analyses. Mon. Not. R. Astron. 
Soc. 2 (2022). Doi: 10.1093/mnras/stac999 

M. Stanbro, M. Briggs, O. Roberts, et al., A study of consec- 
utive terrestrial gamma-ray flashes using the gamma-ray burst 
monitor. J. Geophys. Res.-Space. 123, 9634-9651 (2018). Doi: 
10.1029/2018JA025710 


S. Xiong, Special topic: Gecam gamma-ray all- 
sky monitor. Sci. Sin.-Phys. Mech. Astr. 50, 
129501- (2020). Doi: 10.1360/SSPMA-2020-0457. 


doi:https://doi.org/10.1360/SSPMA-2020-0457 

L. Lin, S. Xiao, Y. Huang, et al., Observational prospects 
for magnetars with gecam. Sci. Sin.-Phys. Mech. Astr. 
50, 129521- (2020). Doi: 10.1360/SSPMA-2019-0397. 
doi:https://doi.org/10.1360/SSPMA-2019-0397 

Y. Huang, Q. Luo, B. Zhang, et al., Ultra-long gamma-ray 
bursts and ultra-soft gamma-ray bursts. Sci. Sin.-Phys. Mech. 
Astr. 50, 129504- (2020). Doi: 10.1360/SSPMA-2019-0415. 
doi:https://doi.org/10.1360/SSPMA-2019-0415 

Y. Su, W. Chen, S. Xiong, et al., Monitoring and research of 
high-energy solar flare emissions with gecam. Sci. Sin.-Phys. 
Mech. Astr. 50, 129505- (2020). Doi: 10.1360/SSPMA-2020- 
0012. doi:https://doi.org/10.1360/SSPMA-2020-0012 


X. Li, X. Wen, Z. An, et al, The gecam and 
its payload. Sci. Sin-Phys. Mech. Astr. 50, 
129508- (2020). Doi: 10.1360/SSPMA-2019-0417. 


doi:https://doi.org/10.1360/SSPMA-2019-0417 

Z. An, X. Sun, D. Zhang, et al., The design and performance of 
grd onboard the gecam satellite. Radiat. Detect. Technol. Meth- 
ods 006 (2022). Doi: 10.1007/s41605-021-00289-y 

D. Zhang, C. Zheng, J. Liu, et al, The perfor- 
mance of sipm-based gamma-ray detector (grd) of 
gecam-c. Nucl. Instrum. Methods Phys. Res. A. 1056, 
168586 (2023). Doi: 10.1016/j.nima.2023.168586. 
doi:https://doi.org/10.1016/j.nima.2023.168586 

C. Zheng, W.X. Peng, X.B. Li, et al., Electron non-linear light 
yield of labr3 detector aboard gecam. Nucl. Instrum. Methods 
Phys. Res. A. . Doi: 10.1016/j.nima.2022.167427 

W. Lu, L. Wang, Y. Yuan, et al., Monte carlo simulation for 
performance evaluation of detector model with a monolithic 
labr3 (ce) crystal and sipm array for ~y radiation imaging. Nucl. 
Sci. Tech. 33, 107 (2022). Doi: 10.1007/s41365-022-01081-3 

C. Zheng, Z. An, W. Peng, et al., Ground calibration of gamma- 
ray detectors of gecam-c. arXiv preprint arXiv:2303.00687 . 
Doi: 10.48550/arXiv.2303.00687 

W. Wang, X. Li, J. Wu, et al., Development and perfor- 
mance study of a dual-layer compton camera (in chinese). 
Nucl. Tech. 46, 030401-030401 (2023). Doi: 10.11889/j.0253- 
3219.2023.hjs.46.030401 

F. Hua, Y. Ye, X. Yang, et al., Several nuclear physics ex- 
periments based on the beijing cyclotron (in chinese). Nucl. 
Tech. 46, 080004—080004 (2023). Doi: 10.11889/j.0253- 
3219.2023.hjs.46.080004 

J. Zhao, S. Jiang, Y. Li, et al., The contribution of 4- 
400 mev cosmic ray protons to the annihilation radiation 
in near-lunar space orbit gamma spectra (in chinese). Nucl. 
Tech. 46, 050201-050201 (2023). Doi: 10.11889/j.0253- 
3219.2023.hjs.46.050201 

M. Moszyski, A. Syntfeld-Kauch, L. Swiderski, et al., 
Energy resolution of scintillation detectors. Nucl. In- 


[24] 


[25] 


[26] 


[27] 


[28] 


[29] 


[30] 


[31] 


[32] 


[33] 


[34] 


[35] 


[36] 


[37] 


[38] 


[39] 


[40] 


[41] 


11 


strum. Methods Phys. Res. A. 805, 25-35 (2016). Doi: 
10.1016/j.nima.2015.07.059 

P. Limkitjaroenporn, W. Hongtong, W. Chaiphaksa, 
et al, The light yield non-proportionality and electron 
energy resolution study of csi(tl) scintillator by comp- 
ton coincidence technique (cct). Mat. Today-Proc. 5, 
15110-15114 (2018). Doi: 10.1016/j.matpr.2018.04.066. 
doi:https://doi.org/10.1016/j.matpr.2018.04.066 

G. Yang, L. Hua, F. Lu, et al., Response functions of a 4 7 
summing ~y detector in G-oslo method. Nucl. Sci. Tech. 33, 68 
(2022). Doi: 10.1007/s41365-022-01058-2 

M. Dong, Z. Yao, Y. Xiao, Development and preliminary re- 
sults of a large-pixel two-layer labr3 compton camera proto- 
type. Nucl. Sci. Tech. 34, 121 (2023). Doi: 10.1007/s41365- 
023-01273-5 

V. Ranga, S. Rawat, S. Sharma, et al., Intrinsic resolu- 
tion of compton electrons in cebr 3 scintillator using com- 
pact cct. IEEE T. Nucl. Sci. 65, 616-620 (2017). Doi: 
10.1109/TNS.2017.2779888 

K. Roemer, G. Pausch, C. Herbach, et al., in IEEE NSS/MIC 
2009, A technique for measuring the energy resolution of 
low-z scintillators. IEEE, 2009, pp. 6-11, doi: 10.1109/NSS- 
MIC.2009.5401909 

L. Swiderski, R. Marcinkowski, M. Szawlowski, et al., Non- 
proportionality of electron response and energy resolution of 
compton electrons in scintillators. IEEE T. Nucl. Sci. 59, 222— 
229 (2012). Doi: 10.1109/TNS.2011.2175407 

S. Zhang, J. Xia, T. Sun, et al., Transition edge sensor-based 
detector: from x-ray to y-ray. Nucl. Sci. Tech. 33, 84 (2022). 
Doi: 10.1007/s41365-022-01071-5 

Hamamatsu (2023). https://www.hamamatsu.com.cn/cn/zh- 
cn/product/optical-sensors/pmt.html. 

X. Han, S. Wang, H. Wu, et al., Csi-bowl: an ancillary detector 
for exit channel selection in y-ray spectroscopy experiments. 
Nucl. Sci. Tech. 34, 133 (2023). Doi: 10.1007/s41365-023- 
01289-x 

Canberra (2023). https://www.mirion.com/products/technologi- 
es/spectroscopy-scientific-analysis. 

BEGe(2023). https://www.gammadata.se/assets/Uploads/BEGe- 
SS-C493 18. pdf 

X. Li, J. Ren, X. Ruan, et al., Particle discrimina- 
tion measurement of liquid scintillators using dt5751. Ann. 
Rep. Chin. Inst. Atom. En. 00, 147-147 (2014). Doi: 
CNKI:SUN: YNXB.0.2014-00-055 

T. Liu, H. Song, Y. Yu, et al., Toward real-time digital pulse 
process algorithms for csi (tl) detector array at external tar- 
get facility in hirfl-csr. Nucl. Sci. Tech. 34, 131 (2023). Doi: 
10.1007/s41365-023-01272-6 

X. Qian, H. Sun, C. Liu, et al., Simulation study on perfor- 
mance optimization of a prototype scintillation detector for the 
grandproto35 experiment. Nucl. Sci. Tech. 32, 51 (2021). Doi: 
10.1007/s41365-021-00882-2 

Y. Wei, M. Guan, W. Xiong, et al., Consistency test of pmt 
spe spectrum from dark-noise pulses and led low-intensity 
light. Radiat. Detect. Technol. Methods 2, 11 (2018). Doi: 
10.1007/s41605-018-0042-6 

M. Li, Z. Wang, C. Liu, et al., Performance of compact plastic 
scintillator strips with wavelength shifting fibers using a pho- 
tomultiplier tube or silicon photomultiplier readout. Nucl. Sci. 
Tech. 34, 31 (2023). Doi: 10.1007/s41365-023-01175-6 

H. Wu, C. Li, A root-based detector test system. Nucl. Sci. 
Tech. 32, 115 (2021). Doi: 10.1007/s41365-021-00952-5 

S. Guo, J. Wu, D. Hou, The development,performances and 


applications of the monochromatic x-rays facilities in(0.218- 
301)kev at nim,china. Nucl. Sci. Tech. 32, 14 (2021). Doi: 
10.1007/s41365-021-00890-2 

[42] D. Hou, J. Wu, S. Guo, et al., The realization and study of 
(21-301) kev monochromatic x-rays. Nucl. Instrum. Methods 
Phys. Res. A. 927,. Doi: 10.1016/j.nima.2019.02.024 

[43] S. Guo, Z. Jiang, J. Wu, et al., Research on a tunable monochro- 
matic x-rays source in (540) kev. Appl. Radiat. Isot. 181, 
110096- (2022). Doi: 10.1016/j.apradiso.2022.110096 

[44] J. He, Z. An, W. Peng, et al., Ground-based calibration and 
characterization of labr3-sipm-based gamma-ray detector on 
gecam satellite: 8-160 kev. Mon. Not. R. Astron. Soc. 525, 
3399-3412 (2023). Doi: 10.1093/mnras/stad2439 

[45] J. Wen, X. Zheng, J. Yu, et al., Compact cubesat gamma-ray 
detector for grid mission. Nucl. Sci. Tech. 32, 99 (2021). Doi: 


12 


10.1007/s41365-021-00937-4 

[46] X. Wen, J. Sun, J. He, et al., Calibration study of 
the gamma-ray monitor onboard the svom satellite. Nucl. 
Instrum. Methods Phys. Res. A. 165301 (2021). Doi: 
10.1016/j.nima.2021.165301 

[47] X. Li, C. Liu, Z. Chang, et al., Ground-based calibration and 
characterization of the he detectors for insight-hxmt. J. High. 
Energy Astrophys. 24,. Doi: 10.1016/j.jheap.2019.09.003 

[48] L. H.R., W. J.J., L. J.C., et al., Lege detector intrinsic efficiency 
calibration for parallel incident photons. Appl. Radiat. Isot. 
109, 551-554 (2016). Doi: 10.1016/j.apradiso.2015.11.102 

[49] M. Alekhin, J. Haas, I. Khodyuk, et al., Improvement of 
c-ray energy resolution of labr3 :ce3+ scintillation detectors 
by sr2+ and ca2+ co-doping. Appl. Phys. Lett. 102,. Doi: 
10.1063/1.4803440 


